Growing Cayley trees described by Fermi distribution 
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We introduce a model for growing Cayley trees with thermal noise. The evolution of these 
hierarchical networks reduces to the Eden model and the invasion percolation model in the limit 
T — > 0, T — > oo respectively. We show that the distribution of the bond strengths (energies) is 
described by the Fermi statistics. We discuss the relation of the present results with the scale-free 
networks described by Bose statistics. 



I. INTRODUCTION 

Recently it has been shown that Bose statistics M can 
be used to describe a scale-free network J|,|| with fit- 
ness of the nodes M . Since scale- free networks are con- 
tinuously growing and develop a power-law connectivity 
distribution it is interesting to investigate their relation 
with self-organized processes l^pQ. In order to address 
this problem we present a model of invasion percolation 
with temperature defined on a Cayley tree, self-organized 
in the limit T — > 0. Our results show that the model 
can be solved analytically using the same technique used 
in the case of the scale-free networks described by Bose 
statistics. 

The invasion percolation model || is the most famous 
and simple example of evolution with quenched disor- 
der. It describes the displacement of a fluid in a porous 
medium. The porous medium is given by a random net- 
work constituted by bonds with different strengths p cho- 
sen with uniform probability in the interval (0,1). The 
classical asymptotic structure generated in this way is a 
fractal and the distribution of the strength values at the 
interface converges in time to a step function 6{p — p c ), 
where p c is the percolation threshold of the static per- 
colation problem. In order to include the effect of fluc- 
tuations on the dynamics of invasion, present in a real 
stochastic cases, we include a temperature-like noise T 
^-|TT| . A structure in which the invasion percolation dy- 
namics can be defined is a Cayley-tree jl2|,[l3| , also known 
as the Bethe lattice. In this structure there are no loops, 
and the number of nodes in the bulk are of the same or- 
der of magnitude as the nodes at the interface. Therefore 
a Cayley tree is considered to be a good representation 
of a d = oo space and it is used in mean field calculations 
[|L4| and in the study of branching processes JTa] . In this 
work we find that the distribution of bond strengths at 
the interface is no more a step function but it is described 
by the Fermi distribution with temperature T where the 
bond strength plays the role of energy. By comparison 
with scale-free networks following the Bose statistics I 
show that both networks grow continuously in time: in 
the power-law network at each time a node is connected 
to the network by m links while in the Cayley tree model 
at each time a node grows giving rise to m new nodes. 
The dynamics of the two networks change in time award- 



ing the fitter nodes in the power-law network or choosing 
the less fit nodes to grow. At the same time in the power- 
law network the distribution of the energies of the chosen 
nodes converges to a Bose distribution while in a Cayley 
tree model the distribution of the energies at the interface 
converges to a Fermi distribution. 
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FIG. 1. Symmetric construction of a power-law network 
and the Cayley tree model considered in this paper. 



II. THE MODEL 

The Cayley tree (or Bethe lattice) is a loop-free net- 
work in which there are three classes of nodes: the root 
node, which is at the origin of the tree and has connec- 
tivity m, the nodes at the interface with connectivity one 
and the nodes in the bulk (below the interface) with con- 
nectivity m + 1. 

We start from the root of the tree (node i = 1) and we 
link it to m new nodes i = 2, 3, . . . , m + 1. We indicate 
each node with a subsequent number, tj indicating the 
time in which it is arrived in the interface. 

At each timestep we choose one node to grow, giving 
rise to m new nodes. Consequently, the interface of the 
tree grows linearly in time, and the growing node is cho- 
sen at each time from the growing number of active ones. 
In order to mimic the quenched noise of the medium we 
assign to each node of the tree an energy e from a fixed 
random distribution p(e). 
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highest energy value are allowed to grow and the model 
reduces to invasion percolation S on a Cayley tree. 




III. EDEN MODEL ON A CAYLEY TREE 
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FIG. 2. Description of the dynamics of the model m — 2. 
At time t = 1 the root node i = 1 with energy ei grows giv- 
ing rise to m new nodes i = 2, 3 with energies £2, £3- Node 
i = 1 is in the bulk or below the interface (connectivity m) 
while the nodes i = 2, 3 have connectivity one and are at the 
interface. At time t — 2 the node i = 2 is chosen to grow, it 
leaves the interface giving rise to m new nodes, i = 4, 5. At 
time i = 3, 4 nodes i = 4 and i = 3 are chosen to grow. In 
the right hand side of the figure we plot the density of states 
of the node in the interface. For example at time t = 4 the 
nodes at the interface are i = 5, 6, 7, 8, 9 and we can draw the 
density of states of the node at the interface by placing a par- 
ticle in the energies levels £5 = £8,£g < £6 < £7 as indicated 
in the picture. 

We assume that higher energy nodes are more likely to 
grow than lower energy ones and that the probability IT; 
for the active node i (with energy ej) to grow at time t 
is given by 



n* = 



E 



jeint(t) 



(1) 



where the sum in the denominator is extended to all 
nodes j that belong to the interface Int(t) at time t. 
The model depends on the parameter /?. Tuning (3 we 
change the nature of the model and the spatial aspect of 
the tree. In the /3 — > limit, high and low energy nodes 
are equally probable to grow and the model reduces to 
the Eden model while in the (3 — > 00 limit the dynam- 
ics becomes extremal such that only the nodes with the 



Let us assume that every node has the same energy 
e°, i. e. p(e) — 6(e — e°). In this case all nodes at the 
interface are equally likely to grow and we call this model 
the Eden model on a Cayley tree. The probability that 
a node i of the interface Int(t) grows at time t is given 
by 



n, 



1 



N Int (ty 



(2) 



where Ni nt (t) is the total number of active nodes. Since 
at each timestep a node of the interface grows, becom- 
ing part of the bulk, and m new active nodes are gener- 
ated, after t timesteps the model generates an interface 
of Nj nt (t) nodes, with 



N Int (t) = (m-l)t + l. 



(3) 



We denote by p(t, ti) the probability that a node, born 
at time ti is still active at time t. Since every node grows 
with probability IT Eq. |^ only if i is a node of the inter- 
face, in mean field p(t,ti) follows 



dp{t,U) _ p(t,ti) 



dt 



N Int (tY 



(4) 



Substituting (g) in (Q) in the limit t — > 00 we get the 
solution 



P(t,ti) 



l/(m-l) 



(5) 



Consequently each node i that arrives at the surface at 
time ti, remains at the surface with a probability that 
decreases in time as a power-law. On the other side the 
same power-law describes also the distribution of the age 
t of the nodes at the interface observed at time t. In 
fact, the probability P(t) that a node born at time r is 
still active at time t, is given by 



P(t) 



rx l/(m-l) 
1. 



(6) 



Thus asymptotically in time the same power-law de- 
scribes the time evolution of the nodes born at time ti, 
(p(t,ti)) and the age distribution of the nodes in the in- 
terface, (P(t)). 
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1. Numerical simulations 

In order to verify the theoretical predictions, we have 
performed numerical simulations of the Eden model on a 
Cayley tree with m = 2,4, 6. In Fig. ||, we report the age 
distribution P(t) of the nodes at the interface for Cay- 
ley trees with m x 10 4 nodes and m — 2,4, 6. The data, 
averaged over 100 runs, follows the power-law predicted 
by (H|). Numerical data are reported together with the 
theoretically predicted power-law (^). 




FIG. 3. Age distribution P(r) of the nodes at the interface 
of a Cayley trees with connectivities m = 2, 4, 6 and 10 4 gen- 
erations. Data have been averaged over 100 runs. The solid 
lines represent the power-laws predicted by (^j) with exponent 
1 - 1/m. 



IV. CAYLEY TREE WITH ENERGIES 

At finite temperature, (3 ^ it is necessary to take into 
account the fact that each node has a different energy 
that defines its dynamics. In this case the probability IT 
for a node i at the interface with energy to grow at 
time t is given by dll) 



n, 



E 



jeint(t) 



(7) 



Since only nodes at the interface can grow, the prob- 
ability that node i would leave the interface at time t is 
given by the product of Pi(t\ei,ti) (the probability that 
the node is active at time t) and IT (the probability that 
the node is chosen to grow in between the active nodes) . 
Consequently, pj(t|ej,tj) decreases in time following 



dp(t\e u ,U) e^i P {t\e ti ,ti) 



at 



Ej 



£7nt(t) 



(8) 



were the sum in the denominator is extended to all the 
nodes j that are part of the interface Int(t) at time t. 

2. Analytic solution 

To solve (|8|) we assume that in the thermodynamic 
limit the sum Z s in the denominator of the l.h.s. of Eq. 
(§), given by 



jelnt(t) 

self averages and converges to its mean value 

J = l,N 



(9) 



(10) 



Moreover, since Z s is an extensive quantity we can self- 
consistently assume 



< z s z F t{i + o{r a )). 



(ii) 



Using (|l l[) , the dynamic equation (|3|) can be written as 

dp(t\e ti ,U) e^p(t\ei,U) 



at 



t 



(12) 



Consequently we found that the time evolution of 
Pi{t\ei,ti) follows a power-law, 



Pi(t\ei,ti) = [ j 



(13) 



but there is multiscaling in the system, i.e. the dynamic 
exponent depends on the energy q of the node, 



m = — 

ZF 



(14) 



After substituting pi(t\ei,ti) from Eq. Jl3| ) with /(e) 
given by Eq. (|l4|), into Eq. (fj"l|), and the sum with an 
integral, we get the self-consistent equation 



<J2^Pj{t\e j ,t j )> '=°° / de p(e) f dt'e 

= z F t{l + 0{t- a ) 

where 

a = (1 + min/(e)) > 1, 



/Be ' 1 \ 



(15) 



zp — m J dep(e) — 
Finally, if we define [ip as 



/3t 



zp = e 



Put 



(16) 
(17) 

(18) 
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we obtain the self-consistent Eq. ( |l7| ) can be interpreted 
as a definition of [i F and it is formally equivalent to 
the definition of the chemical potential in an equilibrium 
Fermi gas 



1 



- = [ dep(e)- 
m J 



■P(z-»f) + 1' 



(19) 



suggesting that many properties of this model can be de- 
scribed by the Fermi statistics. 

On the other hand the probability P(t) that a node 
born at time r is still active at time i, is given by a 
power-law 



P(t)= I dep(e)[- 



(20) 



3. Existence of the solution for the chemical potential 



Equation ([L9|) always has a solution for the chemical 
potential [i F . In fact, since p(e) is a normalized distri- 
bution function, the integral I{p(e), /j, f } on the r.h.s. of 
equation (pf), 



I{p(e),^ F }= / dep(e)n F (e) 



for a fixed fi F is bounded by 
1 

< I{p(e),n F } < 



(21) 



(22) 



implying that for the chemical potential fj, F the solution 
of Eq.(|l9|), I{p(e),[i F } = 1 — 1/m, satisfies 

1 



\ log(ra - 1) < hp < e r 

P 



-log(m-l). (23) 
P 



This proves that for a real tree with m > 1 the equation 
( |l9| ) always has a solution. Finally some attention should 
be given to the special limits (3 — > and (3 — > oo. 

/3 — ► Zimit - In this case we recover the solution of the 
Eden model on the tree, zp — m — 1. Since the probabil- 
ity distribution p(e) is normalizable and the occupation 
number 



n F (e) 



1 



Zp l + l' 



(24) 



equation (|19| ) reduces to zp — > m — 1 and thus /3/ij? — > 
ln(ro — 1) in such a way that /xj? > [zp > 1) if m > 2. 

/? — ► oo limit- In this limit the Fermi-Dirac distribu- 
tion converges to the step function 

n F [e) -» 0(e - mf) (25) 

and the self-consistent equation (|l9|) becomes 



1 

1 = 

m 



(26) 



^. Mass conservation 



The self-consistent relation ( |19[ ) can also be derived 
from mass conservation, i.e. from the knowledge that 
the total number of nodes at the interface is given by 
N ~ (m — l)t. Consequently, 



N=(m-l)t= P(*l £ i>*i)- 

i£[0,mt] 



(27) 



We can substitute the sum in the right hand side of 
Eq. ( p7| ) with the mean over the energies e% of the nodes 
% of generation U. Moreover, in the thermodynamic limit 
we can approximate the sum over i with an integral over 
tii the mass conservation relation becoming 



(m — l)t = m I dep(e) j drp(t\e,T) 
m 



dep(e) I dr (I)' 
mt I dep(e)—7- r ; 



(28) 



where in the last equation we have neglected terms of 
order 0(t~ a ). Thus both the mass conservation relation 
( p7| ) and the self-consistent relation ( |l9| ) allow us to de- 
fine the chemical potential fip , describing the evolution 
of the network as the chemical potential of an equilib- 
rium Fermi gas with specific volume v c = 1 + l/(m — 1). 
However this last expression explains the meaning of that 
relation. In fact, the number Nj nt (e) of nodes with en- 
ergy e at the interface at time t is given by 



N Int (e) = mtn F (e)p(e), 



(29) 



where np(e) is given by the Fermi occupation number 

1 



n F (e) 



e P(n-n F ) _|_ i ' 



(30) 



In other words, the distribution of the energy at the inter- 
face reaches a stationary limit given by ( p9[ ) and defined 
by a Fermi distribution with chemical potential given by 
(|l9|). In the mean time, the density of nodes with energy 
e present in the bulk, NBuik{t), reaches a stationary limit 
as well. In fact, since the nodes in the bulk are the ones 
of the network that are not at the interface, using (|2£ 
we have 



N Bu ik(e) = p(e)[l - n F (e)]. 



(31) 
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5. Asymptotic dynamics 

The dynamical evolution of the network brings the sys- 
tem to the stationary state described by the distri- 
bution function (^0|) , as it has been shown by the solution 
of the dynamical equation (|J). Moreover the dynamics 
stabilizes this distribution. In fact, in the asymptotic 
limit, when the survivability follows dl3| ) the probability 
that a node of energy e will grow and leave the interface 
is given by 

t)=mf de'p(e') f dt' 9 ^^ S(e - e') (32) 



which, using ( |l2] ) can be extimated to be 



iTF(e,t) = / de'mp(e) I dt' 



o 0(e-fj,i 



] p(t\e',t>) 



mp{e)[l - n F (e)]. 



/3(€-M_f) 



S(e-e>) 



(33) 



Consequently, the probability that a node of energy e 
leaves the interface, asymptotically in time, reaches a 
stationary limit independent of the particular evolution 
of the network, given by 



7r F (e, t) TT* F (e) = p(e)[l - n F (e)]. 



(34) 



If we observe an evolving network and we have no knowl- 
edge of the age of the nodes, but only of their energies, 
the complete dynamics is determined by 7r(e, t) describ- 
ing which is the probability that a node with energy e will 
leave the interface at time t. While the complete dynam- 
ics (||) is clearly dependent on time, 7r(e, t), reaches the 
stationary limit n F (e), defining the invariant dynamics 
of the system. 

The stability of the distribution (e) of the ener- 
gies in the bulk, is thus enforced by the dynamics. In fact 
we have found that asymptotically in time the probabil- 
ity that a node with energy e is chosen to grow 7r£,(e) is 
proportional to the number of nodes in the bulk iV^^e) 
given by Eq. (|l|). 



6. Numerical support 



/3 = 2, 5, 10, 20 and compared to the theoretical predic- 
tion Eq.@, P(t) ~ (f ) 5 with S given by e~^ F for the 
uniform distribution p(e) = 1, e G [0, 1]. 




FIG. 4. The energy distribution of the nodes in the inter- 
face for (5 — 5, 10, 30 in the case of a uniform energy distribu- 
tion p(e) — 1 for e € [0, 1], to = 2 and predicted chemical po- 
tential /if = 1/2. The solid lines indicate the predicted Fermi 
distribution. 




Choosing the node energy from a uniform distribution 
p(e) = 1 with e G [0, 1], we have simulated the growth 
of a Cayley tree with m — 2 and various values of (3. In 
Fig. ^ we report the distribution of the energies of the 
active nodes for a network of size N = 2 x 10 4 nodes for 
(3 = 5, 10, 30. The solid line in the figure represents the 
theoretical prediction described by ( p9| ) and ( p0[ ) with a 
chemical potential given by @. In Fig. | the distribu- 
tion of the age of the nodes at the interface is shown for 



FIG. 5. Probability distribution of the age r of the nodes 



at the interface in a Cayley tree with m 
a function of /? = 2, 5, 10, 20. 



2 and time 10 as 



CONNECTION WITH OTHER PROBLEMS IN 
STATISTICAL MECHANICS 
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A. Invasion Percolation 

The static percolation problem on a Cayley tree is 
exactly solved and gives for the percolation threshold 
p c = 1 — 1/m. The invasion percolation problem, rep- 
resents the first example of a self-organized system and 
it has been solved exactly on a Cayley tree in and 
numerically in |l3| . 

If we suppose that the energies of the tree model are 
uniformly distributed between zero and one, i.e. p(e) = 1 
for e E [0, 1], in the limit (3 = oo the evolution of the in- 
terface of the Cayley tree network maps exactly into the 
front of invasion percolation on a Cayley tree. We can 
find the chemical potential /j,f by solving 



1 — l/m= J denp(e) 



ck- 



1 



o e 

de = /zf 

o 



(35) 



Consequently, asymptotically in time, the density of 
nodes with energy e found at the interface follows the 
step function 



n F (e)=0(e-(l-l/m)), 



(36) 



predicting the correct threshold p c = 1 — 1 jm for invasion 
percolation on a Cayley tree. 



B. Percolation 

The percolation transition can be investigated choos- 
ing a distribution function 



p(e) = pS(e-l) + (l-p)6(e), 



(37) 



with A = (p — p c ) and p c = 1/m. But only one solution 
is possible, because the fugacity is strictly positive. This 
implies that if A is positive (p > p c ) we should take the 
negative sign in equation (|3^) while we have to use the 
positive sign for p < p c . 
Thus we have, 



Zpe 



-13 



e^(l-p c )/|A| p<p, 
|A|/p c p > p, 



Consequently, using (|3 
at the interface is: 



(40) 

the mean energy of the nodes 



1 



1 



de e p(e)njr(e) 
V 



1 



Pc e^z'r 



1 



(41) 



In the T — + limit e mean displays a sharp transition at 
p c = 1/m. Using the solution (|40|) in we obtain 



f e ^ppc/IAI p < p c 
\ A/(l-p c ) p>p c 



(42) 



lo- 



ur 



10 
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describing the fact that each node of the Cayley tree has 
energy e = 1 with probability p and energy e = with 
probability 1 — p. In the limit T — > only the nodes 
with e — 1 would grow. Consequently, assigning to the 
energy at the interface a meaning of cost function 
we have that in the T = the mean energy of the node at 
the interface is minimized. When we consider the proba- 
bility distribution (|37|), the self-consistent equation jl9| ) 
reduces to a quadratic equation for the fugacity zp, 



P 



°-Pz~ F x 



that in the (3 — > oo limit has the solutions 



-A 



2(1 -Pc) 



1± 



l + 4e -^ c(1 - pc) 



A 2 



(38) 



(39) 



FIG. 6. Mean energy of the nodes at the interface as a 
function of p and j3 = 10, 5, 3, 1. The solid line represent the 
theoretical prediction in the j3 — > oo limit. 

In order to verify the predicted solution we have sim- 
ulated the growth of a Cayley tree with m = 3 branches. 
We have assigned to each node an energy equal to or 
1 with probability given by fl37|). The evolution of the 
network described by (^) proceeds preferentially on high 
energy nodes, in such a way that in the (3 — ► oo the mean 
energy of the interface is minimized. In Fig. |(j we report 
the mean energy of the nodes at the interface as a func- 
tion of p for different temperatures. As the temperature 
goes to zero, the mean energy of the system has a sharp 
transition from zero to a non zero mean energy distribu- 
tion. This shows that during the network evolution the 
nodes with high energy remain at the interface. 
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VI. PARTITION IN TIME, ENERGY AND 
MEASURE 



A. Partition in time 

If we coarse grain the time and divide the nodes fol- 
lowing their age, we observe that the nodes with approx- 
imately the same age have a survivability p(t,U) that 
grows in time as a power-law of ii = t/ti, 



< p(t,ti) > e 



dep(e)p(t,U 



it 



with Tp asymptotically given by 



(43) 



(44) 



This power-law describes the fact that the probability to 
be at the interface is lower for older nodes, but decreases 
much faster for nodes of energy e > pp. The scale of the 
power-law depends on the ratio between the observation 
time and the time U when the node i joins the interface: 



ti = -■ 

U 



(45) 



Because of this structure p(t,ti) is not a function sepa- 
rately of t and U, and not even of their difference t — ti, 
but only of their ratio i^. Thus ( (43| ) describes both the 
probability that a node born at time ti is at the interface 
at time t and the probability that a node at the interface 
at time t is arrived in the network at time ti. 

Consequently the network at a given time will contain 
nodes of every age at different stages of their evolution. 
In other words the tree acts as if it was recording the evo- 
lution of the nodes, and the time partitions of the Cayley 
tree behaves as riscaled copies of the same network. 



B. Partition in energy 

A different case implies that the ages of the nodes 
are unknown but we have information about their en- 
ergies. Thus it is natural to group together nodes with 
the same energy. 

Since the total number of nodes with energy e in the 
Cayley tree are p(e)t, the fraction of the nodes that re- 
main at the interface is given by the integral of p(t\U, e) 
over ti, normalized to p(e)t, i.e. 



1 



p(e)t 



dtip(e)p(t\ti,e) = n F (e). 



(46) 



As observed above this density is given by the Fermi oc- 
cupation number np(e), reached in the thermodynamic 
limit, and it gives the characteristic stable distribution 
of the Cayley tree evolution. 



C. The partition functions 

The emergence of the quantum statistics in the de- 
scription of the geometrical structure of the interface of 
the Cayley tree suggests us the possibility to perform a 
more detailed investigation of its statistical properties. 



1. Static partition function 

In order to study the aspect of the Cayley tree at a 
given time, we can define the static partition function 
Z s (t) as 



Z s (t) 



E 



ie/nt(t) 



(m - 1) 



(47) 



This describes the distribution of the energy values in 
the interface. As proven previously in (0), the static 
partition function, asymptotically in time, reaches a sta- 
tionary limit and satisfies 



lim 



z s {t) 



z F 



(m — 1) 



(48) 



Consequently, the static partition function defined in 
Eq. ([!?]) is an extensive quantity related to the chemi- 
cal potential of the network by 



1 lim toEjZSjt)) 

(3 t^oo t 



pp-T\ag(m-\). 



(49) 



2. Dynamic partition function 

In order to describe the Cayley tree evolution we asso- 
ciate to it the sequence {e^} of the energy of the nodes 
selected to grow at time t = r. If we consider all the time 
dependent sequences associated to different realization of 
the networks for t timestcps, we can associate with them 
a dynamical partition function Z D (t) given by 



Z D (t)= J nr=i,t{der}e /J ^=i e -P(ei,e 2 , 



■ e t ) (50) 



where the sum goes over all the different realizations of 
the network and -P(ei, £2, £mt) is the probability of the 
sequence {e,.}. As we have shown previously (Eq. |34|), in 
the thermodynamic limit the probability that an energy 
e will be selected to grow reaches a stationary value given 
by the quantum occupation numbers. Consequently the 
partition function, asymptotically in time can be written 
as 



Z£(f + l)~Z£(t) / dee^TT F (e) 



(51) 
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giving in the thermodynamic limit, 

z D {t) - (z9Y 

were 

z$ =m [ dep{e)e 0e [l - n F {e)]. 



(52) 



(53) 



Thus for the induced dynamic in the energy space we 
can introduce a dynamical partition function describing 
the possible realization of the networks. The dynamical 
partition function, asymptotically in time is given by 



lim 

t — >oc 



log(Zj?(*)) 



t 



iog(*£) 



(54) 



The connection between the dynamical and the static 
fugacity is easily shown to be, 



(m- l)z F + z\ 



to < e p > 



(55) 



We have solved analytically the model studying its char- 
acter in particular on the limit (3 = and then at finite 
temperature. The distribution of strength bonds follows 
a Fermi distribution and the dynamics asymptotically in 
time replicates and stabilizes this distribution. The bond 
strength plays the role of energies in the Fermi distribu- 
tion. The distribution of ages of the node at the interface 
follows an effective power-law. We have then studied the 
cases in which the model reduces to the study of Invasion 
Percolation and percolation. We have found the perco- 
lation threshold and the mean energy of the growth. Fi- 
nally we investigate the statistical properties of the solu- 
tion and we define two partition functions describing the 
ensemble of growing Cayley trees. Finally this system is 
a symmetric construction of a power-law network follow- 
ing Bose distribution, as shown in Fig.l and it opens the 
way to understand the self-organized nature of scale-free 
networks. |fl9| . 



3. Entropy 

One characteristic feature of this model is that the dy- 
namics introduces a distinction between different nodes 
of the tree. In fact the nodes of the tree are distinguished 
between nodes in the interface and nodes in the bulk. The 
separation of this 'two phases' is forced by the dynamics 
itself and can be the reason for the statistical mechanics 
properties described by quantum statistics. In fact, if we 
look at the structure of the fermionic entropy S , 



S 



dep(e)n F (e) \og(n F (e)) + 



dep(e)[l-n F (e)]log([l-n F (e)]) 



(56) 



we can observe that it can be interpreted as the sum of 
two Shannon entropies of the tree: the bulk and of the 
interface entropy, 



S F — Si, 



Sbulk- 



(57) 



VII. CONCLUSIONS 

In this work we have introduced a model for a growing 
Cayley tree with thermal noise characterized by 

• Growth: At each time exactly to nodes are added 
and one is eliminated at the interface, the number 
of nodes in which percolation can occur grows in 
linearly in time as N = (to — l)t nodes; 

• Time dependent dynamics: Each node can perco- 
late only once and the probability for a node to be 
chosen as the percolating one is a decreasing func- 
tion of time. 
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